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1. Introduction 

It is believed that QCD undergoes phase transitions as a temperature T and a 
chemical potential y, increase. For example, an expected phase structure is found 
in Fig. 1 of ref. To theoretically investigate what really happens at non-zero T 
and n in QCD, lattice calculations seem ideal. As reviewed in ref. lattice QCD 
investigations are indeed very successful at finite T but fi — 0. On the other hand, 
lattice QCD simulations at non-zero /i have remained to be extremely difficult 
for a long time, since a complex nature of the lattice QCD action makes a direct 
application of straightforward Monte-Carlo methods impossible. Recently, however, 
several new techniques have been proposed to overcome this difficulty. In this talk 
I mainly review recent progresses in lattice QCD at non-zero it. 

2. Lattice QCD at non-zero /x 

A possible solution to the complex action problem at non-zero it was proposed in 
ref. 03 where a reweighting method from \i = was employed. This method, called 
the Glasgow method, however failed since an overlap between gauge ensembles at 
/i = and at /i ^ is exponentially small in the volume. 

Recently two new methods have been proposed for small [i near the critical 
temperature T c . One is a variant of the Glasgow method 3, : Instead of reweighting 
only /i from /i = 0, a multi-parameter reweighting in both /i and T from it = and 
T = T c (p = 0) is cmployed ^^ ' ' ™ The other method has empl oyed an imaginary 
chemical potential /i = ifii to avoid the complex action problem "^ " I Calculations 
in both methods so far are restricted to N t = 4, where the temperature is given by 

r = i/(jv t o). 

a For a reweighting of only fi by the Taylor expansion, see ref. HI 
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2.1. Multi-parameter reweighting method 

I first discuss the multi-parameter reweighting method, in which the QCD partition 
function at non-zero fi is rewritten as 



where U represents a gauge configuration, j3 = 6/g 2 , m is a quark mass and /i 
is a (quark) chemical potential. Here Sq is a gauge action and det D is a quark 
determinant. In the actual calculation, gauge configurations are generated at fi' = 
and (3' ~ /3 c (m, 0), where (3 c {m, /i) is a critical coupling at /i and to. Since both (3 and 
are used in the reweighting, this method is called a multi-parameter reweighting. 
Once configurations are generated at a' and /3', one can treat the term in {■ • ■} as 
a part of observables. 

Fodor and Kafir*- 1 have employed the multi-parameter reweighting method start- 
ing from gauge configurations generated by 2+1 flavor KS quark action with the 
plaquette gauge action at u' = and f3' = (3 C . Since they have exactly evaluated 
det D(U, fi, to) appeared in the reweighting factor, lattice volumes in their calcu- 
lations are restricted to x 4 with L s = 6, 8, 10, 12. Dynamical quark masses 
in the simulation are m u( i/T = 0.0368 (m^/nip — 0.188(1)) and m s /T = 1.0 
(m^/rriK — 0.267(1)), which are very close to experimental values, m^/nip = 0.179 
and m^/nix = 0.267. 

In order to distinguish the first order phase tran sitio n from the crossover, they 
investigated a volume dependence of Lee- Yang zero /?o, which are zeros of the 
partition function Z{j3, /j,, to) in the complex (3 plane. Since there is no real phase 
transition in the finite volume, Im/?o ^ at V = ^ oo. If limv/^oo Im Pq = 0, 
there exists the first order phase transition at (3 — (3o, while if limy^oo Im/3o =/= 0, 
only a crossover appears at (3 — Re/3o- In the left of Fig^] taken from ref.[71 Im/3^°, 
obtained by extrapolating (3$ to the infinite volume as (3q(V) — (3^ +c/V, is shown 
as a function of fia. The critical endpoint, which separate the first order phase 
transition line from the crossover line, is estimated to be /ia ~ 0.18. Converting 
(Re /3(j° , /ia) to the physical unit (T c ,/^s = 3/i), where /ib is a baryonic chemical 
potential, a phase diagram for QCD at non-zero T and /i is given in the right 
of Fig^ also taken from ref. and the critical endpoint becomes (Te,^%) = 
(162(2), 360(40)) Me\d. 

Instead of evaluating det D(U, /i, to) exactly, Bielefeld-Swansea group has calcu- 
lated it by the Taylor expansion, 




In 



det D(U, /j,, to) 



RiH + i?2M 2 H 



(3) 



det D(U, 0, to) 
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Fig. 1. (Left) Im /3q° a s a function of the chemical potential. (Right) The phase diagram in 
physical unit. 



to perform the simulation on a relatively large volume, 16 3 x 4. Gauge configurations 
are generated by the improved gauge action and the Nf — 2 p4-improved KS 
quark action at m/T — 0.4, 0.8, corresponded to m^/rrip ~ 0.7, 0.85. The transition 
temperature was estimated from the peak position of susceptibilities x^pip an d XL, 
where L is the Polyakov loop. The result of the critical temperature T c (/j,) as function 
of the quark chemical potential /j, q can be found in Fig. 16 of ref. [H] Errors of the 
reweighting remain small at ^7 < 0.4, but they increase as /i q becomes larger. 



2.2. Imaginary chemical potential 

Since the quark determinant det(U, ifii, m) becomes real for the imaginary chemical 
potential ijjij, direct Monte-Carlo simulations are possible in this cas t) I I One then 
fits the critical temperature obtained at the imaginary chemical potential in the form 
of the Taylor expansion: 

T c (im)=T c (0) + c 2 ^ + O(fij), (4) 

as seen in Fig. 9 of ref. ^1 for Nf = 2 and in Fig. 6 of ref. for Nf = 4. Assuming 
an analyticity at p, = 0, one immediately obtain the critical temperature at real /i: 

Tjji) =T c (0)~c 2f i 2 + O(n 4 ). (5) 

In this method, a restriction that |/2/ = ^-| < § exists, due to periodicity and 
symmetry, which leads to 

Z(fii = 3- + x ) = z iPi = 3 -x), (6) 
as explicitly confirmed in Fig. 6 of ref. 1101 
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Table 1. Values of C from various simulations. Here Imag, MulPar- 
Rew, MulParRewTay represent the imaginary chemical potential, the 
multi-parameter reweighting and the multi-parameter reweighting with the 
Taylor expansion, respectively. 



N f 


C 


m ud /T m s /T 


quark 


method ref ■ 


2 


0.0056(4) 
0.008(3) 


0.1 

0.4,0.8 


KS 
p4 


ImagED 

MulParRewTajlHnil 


2+1 


0.0032 


0.0368 1.0 


KS 


MulParRewEI 


3 


0.0028(7) 

0.013(5) 

0.00678(10) 


0.1 
0.02 
0.1-0.16 


p4 
p4 
KS 


MulParRewTaySH 
MulParRewTajCH 
ImagH^ 4 included) 


1 


0.011 


0.2 


KS 


Imagl 



2.3. Comparison of various results 

The critical temperature T c can be parameterized as a function of the baryonic 
chemical potential /is as 

T c{Hb) _ r (J? B 

r c (o) t c (o) 2 ' [) 

Results for C from various simulations are accumulated in tabled which indicates 
that C becomes larger as Nf increases. Further investigations including the quark 
mass dependence, however, will be required for the definite conclusion on the Nf 
dependence of C. 

A critical endpoint (h^,Te) is a point on the critical temperature line which 
separates the cross-over from the first order phase transition. Results of the critical 
endpoint from two groups are given in table [21 from which one may notice that the 
location of the critical endpoint strongly depends on the quark mass. Again further 
investigations will be required for the reliable estimate on the location of the critical 
endpoint for the physical case. 

One may convert the mass dependence of the critical endpoint to the fi depen- 
dence of the critical quark mass as 

mM_ 1 + D fj^\\ (8) 



m c (0) \7rT ( 

where the critical quark mass separates the first order phase transition from the 
cross-over: the phase transition is the first order at to < m c while it becomes cross- 
over at m > m c . For the Nf = 3 case, The imaginary chemical potential methocP^, 
using the KS quark action, gives D = 0.84(36) and m c (Q)/T — 0.123(1), while the 
multi-parameter reweighting method with the Taylor expansion^! using the p4- 
improved KS quark action, leads to D = 690 and to c (0)/T = 0.0028(16). Although 
errors are large and quark actions are different, a discrepancy between two results 
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Table 2. The critical endpoints (fj,g,TE), together with the critical temperature T c at fig = 0. 



N f 


Atf MeV 


T E MeV 


T c MeV 


m ud /T 


m s /T 


quark 


mcthod rof ■ 


2+1 


725(35) 
360(40) 


160(4) 
162(2) 


172(3) 
164(2) 


0.1 

0.0368 


0.8 
1.0 


KS 
KS 


Lee- Yang zer<0 
Lee- Yang zercd 


3 




156(30) 






0.0028(16) 
0.02 


0.0028(16) 
0.02 


p4 
p4 


Binder cumulantj^l 
Binder cumulantii^J 



is huge. We will have to understand and resolve this discrepancy, in order to make 
the definite conclusion for the phase structure at non-zero chemical potential. 



3. Some remarks 

3.1. Reweighting vs. imaginary chemical potential 

One can check a reliability of the reweighting method by comparing it with the 
direct calculation in the case of the imaginary chemical potential. In the left of 
Fig. |21 taken from ref. |SJ {tptp) is plotted as a function of /i/. The multi-parameter 
reweighting method agrees with the direct calculation, while the Glasgow method 
does not. This comparison may suggest that the multi-parameter reweighting is 
more reliable. As mentioned in the previous subsection, however, the result must be 
symmetric in /i/ at fij = n/12 ~ 0.26. This symmetric property is clearly violated 
in both mult-parameter reweighting method and the direct calculation. It is likely 
that the violation may be caused by the insufficient statistics for both methods. 
Therefore, in the case of the reweighting method, one should always check whether 
the symmetric property is satisfied or not for the imaginary chemical potential, 
using the same ensemble employed for the reweighting of real /i. 
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Fig. 2. (Left) (iptp) as a function of Im fj, for the direct simulation(squares), the multi-parameter 
reweighting (crosses) and the Glasgow method(circles). (Right) Im f3a as a function of 1/V at 
fj,a = 0.16, together with the linear fit(solid line), the quadratic fit(dotted line) and the linear fit 
for larger 3 volumes(dashed line). 
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3.2. Cautions to the Lee- Yang zero analysis 

It has been pointed out that the partition function becomes zero in the infinite 
volume limit for all /i ^ 0: lim ^ 0, V) = 0, due to the sign problem of 

V— >oc 

the complex actional. Therefore one should check the volume dependence of not 
only the 1st Lee- Yang zero but also the 2nd, 3rd, • ■ • Lee- Yang zeros, in order to 
distinguish the first order phase transition from the cross-over. 

[71 

In addition to the above subtlety, the procedure of taking V — > oo limit may 
cause some systematic uncertainties in the estimation of the critical endpoint. For 
example, an extrapolation linear in l/V gives Im (3^ = 0.0018(7) at /ia = 0.16, 
suggesting that this point corresponds to the cross-over. As shown in the right of 
Fig. [21 however, the quadratic fit or the linear fit for larger 3 volumes gives Im 
/3^° = —0.0004(11) or 0.0003(9), respectively, concluding that the phase transition 
is consistent with the first order. Clearly a more careful finite volume analysis is 
required for a definite conclusion on the location of the critical endpoint. 

4. Conclusions 

By recent developments for lattice QCD techniques, numerical simulations become 
possible at small ijl and non-zero T. Preliminary results suggest ^ ~ 350 - 450 MeV 
at physical light and strange quark masses. However further confirmations will be 
definitely required for the reliable estimate of /i^. So far numerical simulations at 
non-zero /x have been performed only with KS fermions. Therefore new calcula- 
tions by other fermion formulations such as Wilson/clover or domain- wall/overlap 
fermions will be necessary to check the present results by KS fermions. Definitely 
new ideas will be needed to explore QCD phase structure at large fi and low T. 
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